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Abstract 



We describe the harmonic functions for certain isohedral triakis soUds. They are the 
first examples for which polyhedral harmonics is strictly larger than group harmonics. 
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1 Introduction 



A harmonic function on R" is a distribution solution to the Laplace equation 



Af = {dl + --- + dl)f = 0, 



di = d/dx, 



(1) 



which is necessarily a smooth function. A classical theorem of Gauss and Koebe states that 
a continuous function on M" is harmonic if and only if it satisfies the mean value property 
with respect to the sphere S^~^. On the other hand the Laplacian A admits an invariant- 
theoretic interpretation relative to the orthogonal group 0{n). Namely the M-algebra of 0{n)- 
invariant polynomials in x = (xi, . . . , Xn) is generated by the squared distance function ^p{x) = 
xl + ■ ■ ■ + x'^ and the Laplacian A = ip{d) is the result of substituting d = {di, . . . , dn) into 
^{x). These two characterizations of harmonic functions allow us to generalize the notion in 
two directions, that is, to polyhedral harmonics and to group harmonics. 

Let P be an n-dimensional polytope in M" and P{k) the /c-dimensional skeleton of P, that 
is, the union of all fc-dimensional faces of P. A continuous function f{x) on MP' is said to be 
P{k)-harmonic if it satisfies the mean value property with respect to P{k), that is, if 



where /x^ is the fc- dimensional Euclidean measure on P{k) with |-P(/c)| := ^k{P{k)) being its 
total volume. Let 'Hp[k) denote the set of all P( A;) -harmonic functions on M". A general result of 
Iwasaki p], Theorem 1.1] states that %p{k) is a finite-dimensional linear space of polynomials. 
It is naturally an ]R[(9]-module, that is, stable under partial differentiations. If G{k) is the 
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symmetry group of P{k), then it is also an R[G'(/c)]-module, that is, stable under the natural 
action of G{k). Moreover, if P{k) enjoys ample symmetry, that is, if G{k) acts on M" irreducibly 
then 1-Lp[k) is a finite-dimensional linear space of harmonic polynomials O Theorem 2.4]. 

Let G be a subgroup of 0{n) and M[a;]*^ the M-algebra of G-invariant polynomials of a; = 
(xi, . . . , Xn)- A G-harmonic function on M" is a distribution solution to the system of PDEs: 

ip{d)f = (V^)eMN^), (2) 

where M[x]^ is the maximal ideal of ]R[a;]'^ consisting of those (p{x)^s without constant term: 
(p{0) = (see Helgason [3| Chap. HI]). Let Tic denote the set of all G-harmonic functions on 
M". To define system ([2]), the polynomial if{x) may not range over all IR[a;]^, but only over a 
set of generators of M[x]^. For example, if G is the entire 0{n) then ]R[x]^j?^"'' is generated by 
the squared distance function only, so that system ([2]) reduces to the single Laplace equation 
dl]). In this article, G will be a finite group in general and a finite reflection group in particular. 
Steinberg |i9j shows that if G is a finite group then Tic is a finite-dimensional linear space of 
polynomials with dim'Hc > |G|, the order of G. When G is a finite reflection group, he goes on 
to determine T-Lg explicitly. Namely, as an M[9]-module T-Lq is generated by the fundamental 
alternating polynomial Ag{x) of G, and as an M[G]-module it is the regular representation of 
G, in particular dim^G = |G|. 

When G = G{k) is the symmetry group of P{k), it makes sense to compare 'Hp(k) with 
T-Lcik)- According to a result in [S^, formula (2.14)] there is always the inclusion 

naik) C 'Hp(k) {k = 0,...,n). (3) 

It is known that 'Hp{k) coincides with l-iGik) when P is any regular convex polytope with center 
at the origin, in which case G{k) is an irreducible finite reflection group so that T~Lp(^k) — T~^G{k) is 
determined by Steinberg's theorem; see Iwasaki [S], Theorem 4.4] and the references therein. One 
can show that the coincidence also occurs, for example, for the truncated icosahedron [5, 6, 6] 
of an Archimedean solid, or the soccer ball. So far, however, no polytope with ample symmetry 
has been known for which inclusion (|3]) is strict (for some A; > 0), that is, polyhedral harmonics 
is strictly larger than group harmonics. The aim of this article is to present the first examples 
of such polytopes. They arise from a one-parameter family of isohedral triakis tetrahedra as 
well as from a one-parameter family of isohedral triakis octahedra in three dimensions. The 
former family contains a desired example for = 1, but none for = 0, 2, 3 (see Theorem 12. ip . 
while the latter family contains such an example for every A; = 0, 1, 2, 3 (see Theorem 15. ip . In 
this respect the latter family is more interesting than the former, but in any case we begin with 
the simpler case of the former family and then proceed to the latter. 

2 A Family of Isohedral Triakis Tetrahedra 

An isohedral triakis tetrahedron is obtained from a regular tetrahedron by adjoining to each 
face of it a pyramid (based on that face) of appropriate height, or excavating such a pyramid, 
where a polyhedron is said to be isohedral if its symmetry group is transitive on its faces, that 
is, if all faces are equivalent under symmetries of the polyhedron; see Griinbaum and Shephard 
[3]. Our triakis tetrahedron P depends on a positive parameter r. To describe it more neatly, 
let T be the regular tetrahedron with which we get started; it is centered at the origin O and 
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Figure 1: Adjoining a pyramid to a face of a regular tetrahedron; Od : Od' = r : 1. 

with four vertices A, B, C, D. If d' denotes the center of the face ABC, then the pyramid based 
on this face has the top vertex d that hes on the ray emanating from O and passing through 
d' in such a manner that the distance ratio Od : Od' is r : 1 (see Figured]). The polyhedron P 
has four six-valent vertices A, B, C, D and four three-valent vertices a, b, c, d, where a, b, c 
are defined in a similar manner as the vertex d is, and P has twelve faces and eighteen edges. 

Let us look more closely at the polyhedron P for various values of r (see Figure [2]). The 
values r = 1 and r = 3 are special in that as a point set, P degenerates to the tetrahedron T at 
r = 1 and it becomes a cube C at r = 3. For r > 1, P is obtained from T by adjoining a pyramid 
to each face of T, while for < r < 1, it is obtained by excavating such a pyramid. Moreover 
P is convex if and only if 1 < r < 3, in which interval the value r = 9/5 is distinguished in that 
P becomes a Catalan solid [H [Z], or an Archimedean dual solid, whose dual is the truncated 
tetrahedron [3,6,6]. When r = 1, although the four points A, B, C, d are coplanar, we think 
of ABd, BCd, ACd as distinct faces of P, and thus Ad, Bd, Cd as edges, and d as a vertex 
of P. Similarly, when r = 3, although the four points A, B, c, d are coplanar, we think of 
ABc and ABd as distinct faces of P, and thus AB as an edge of P. With this convention, as a 
combinatorial polyhedron, P has the constant skeletal structure for all r > 0. 

Let G{k) be the symmetry group of P{k). When k = 1, for every r > the group G{k) 
is always the same as the symmetry group of the tetrahedron T, which we denote by WlA^) 
since it is a Weyl group of type A3. On the other hand, when k = 0,2, 3, the group G{k) stays 
the same as W^A^) as long as r 7^ 3, but at r = 3 it jumps up to be the symmetry group of 
the cube C, which is denoted by iy(B3) as being a Weyl group of type B3. For simplicity of 
notation, TiwiAa) and T-LwiBs) are abbreviated to and T-Lb^ respectively. Note that 

"^Ag £ "^Bs, dimT^Aa = 24, dim^Bg = 48. 

It is our problem to consider how T-Lp^k) behaves as r varies and when l-Lpi^k) is strictly larger 
than T-Lcik)- Any value of r with this phenomenon is referred to as a critical value. It turns out 
that the vertex {k = 0) and face {k = 2) problems have no critical values, but the edge {k = 1) 
problem certainly has a critical value ri = 3.62398 ■ ■ ■ . It is an algebraic integer of degree six 
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0<r<l 0<r<l 

Figure 2: A one-parameter family of isohedral triakis tetrahedra. Each polyhedron is shown 
in a "cardboard" view on the left, as well as in a skeletal view on the right. Polyhedron with 
1 < r < 3 is not indicated, but imaginable from the one with r > 3, although the former is 
convex while the latter is not. See also a series of pictures in Griinbaum p], Fig. 14]. 
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that is the unique positive root of a sextic equation 

Xi{r) := + 2r^ + - 36r^ - ibr^ - 270r - 405 = 0. (4) 

How this algebraic equations arises will be explained in §4.21 The volume problem {k = 3) need 
not be dealt with independently, since the isohedrality of P implies 'Hp(2) = '^p(3) by a result 
of [SI Theorem 2.2] so that the face and volume problems have the same solution. 

Theorem 2.1 For the family of isohedml triakis tetrahedra, 1-Lp(k) 'is strictly larger than T-Lcik) 
if and only if k = 1 and r = ri = 3.62398 ■ ■ ■ . The function space T-Lpi^k) 'is given by 

[Hbs {if either k = 1, r = ri; or k = 0,2,3, r = 3), . . 

I rtAg [otherwise) . 

This theorem asserts that when k = 1 and r is at the critical value ri, the figure -P(l) has 
only the same symmetries as a tetrahedron, but 'Hp(i) jumps up to become the space of cubic 
harmonics, which is strictly larger than that of tetrahedral harmonics. For each = 0, 2, 3, a 
jumping phenomenon also occurs at r = 3 with the space 'Hp{k) jumping from "Has to "Hbs, but 
at the same time the group G{k) also jumps from W^A^) to 1^(63). Altogether, the equality 
^p(fc) = 'HG(fc) continues to hold and no critical phenomenon occurs at r = 3. After a review 
in ^ on PDEs that characterize polyhedral harmonics. Theorem 12.11 will be established in §H 

3 Invariant Differential Equations 

Iwasaki [SI Theorem 2.1] derives a system of partial differential equations 

T^\d)f = (m = 1,2,3,...) (6) 

that characterizes 'Hp(k) as its distribution solution space. Here is a brief review of it when P 
is a three-dimensional polyhedron and = 0, 1, 2. For j = 0, 1, 2, let {Pj^ jj^g/^ be the set of all 
j-dimensional faces of P, where Ij is an index set. Let Hi be the j-dimensional affine subspace 
of M? containing Pi.. Let Pi. be the foot of orthogonal projection from the origin O down to 
Hi-. We mean by ij -< ij^i that Pj^. is a face of Pi^+i- If ij -< ii^^^ then the vector pi,. — pi_.^^ is 
parallel to the outer unit normal vector of dPi in Hi at the face Pj , so that a number 

[ij : ij+i], called the incidence number, is defined by the relation pi. — Pi^^^ = [ij : ij+i]nj^,jj+i. 
Put /(I) = {i = {io,ii) '■ io -< ii}, /(2) = {i = {io,ii,i2) '■ «o -< ^1 -< ^2}, and define 

[io'-ii] {i = {io,ii) e 1(1)), 

[io : ii][ii : ^2] (i = (^0,^1,^2) e 1(2)). 



Then Tm\x), k = 0,1, 2, are homogeneous polynomials of degree m defined by 

r^\x):=J2{p,„xr, (7) 

-r^K^) ■= Yl ^]hrn{{Pio,x), {Pi„x)), (8) 

ie/(i) 

Tmi^) ■= i^^rn{{Pio,x), {pi„ x) , {pi^,x)), (9) 

j6/(2) 
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Figure 3: Combinatorial data for the system ([6]) of PDEs. 



where ( ■ , ■ ) is the standard inner product on M.^ and hm stands for the complete symmetric 
polynomial of degree m in two or three variables (see e.g. Macdonald [8J). 

The general construction mentioned above will be applied to the particular case of our 
polyhedron P upon adjusting the notation to the current situation. In order to represent the 
index sets Jq, h, h, it would be best to let the vertices, edges and faces to speak of themselves: 

Jo = {A, B, C, D} U {a, b, c, d}, 

Ji = {AB, AC, AD, BC, BD, CD} U {Ab, Ac, Ad, Ba, Be, Bd, Ca, Cb, Cd, Da, Db, Dc}, 
I2 = {ABc, ABd, ACb, ACd, ADb, ADc, BCa, BCd, BDa, BDc, CDa, CDb}, 

where {■ ■ ■ } stands for an orbit under symmetries. For an index Ad G /i, the same symbol 
Ad denotes the foot of orthogonal projection from the origin O to the affine line ^Ad passing 
through A and d; this rule also applies to another index AB G /i as well as to all the other 
indices of Ji. In a similar manner, for an index ABd G I2, the same symbol ABd denotes the 
foot of orthogonal projection from O to the affine plane -ffABd passing through A, B, d, with 
this rule applying to all the other indices of I2; see Figured 

Up to symmetries there are three types of vertex-to-edge incidence relations: 

(VEl) A^AB, (VE2) A ^ Ad, (VE3) d ^ Ad. (10) 

On the other hand, up to symmetries there are two types of edge-to-face incidence relations: 

(EFl) AB^ABd, (EF2) Ad ^ ABd. (11) 

The incidence number of an incidence relation depends only on its type. Let ve,^ denote the 
incidence number of type (VEz/), u = 1,2,3, and efj, denote that of type (EFu), u = 1,2. 
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For example we have [A : Ad] = ve2 and [Ad : ABd] = ef2, so that A — Ad = ve2 ?^A,Ad and 
Ad — ABd = ef2 riAd.ABd for A -< Ad -< ABd. These numbers are positive in the situation of 
Figure [3|, but they or some other incidence numbers may be negative for some values of r. 
For explicit calculations, it is convenient to work with Cartesian coordinates by putting 

A = (1,-1,-1), B = (-1,1,-1), C = (-1,-1,1), D = (1,1,1). (12) 

The symmetry group W^(A3) of the tetrahedron T then admits an invariant basis 

/ \ 2i2i2 / \ / \ 22i22i22 

62(0;) := + X2 + X3, 63(0;) := X1X2XS, e4[x) := X2X^ + x^x-^ + X1X2, 

so that T-La-, can be characterized as the solution space to the system of PDEs: 

e2id)f = esid)f = e,id)f = 0. (13) 

As an M [(9] -module, T-La^ is generated by the fundamental alternating polynomial 

Z\a3(x) = {xl - xl){xl - xl){xl - xj). 

If the vertices of T are taken as in formula f[T^ . then C is the cube with vertices at 
(±1,±1,±1). The symmetry group iy(B3) of C then admits an invariant basis 62(0;), e4(x), 
66 (x) := e|(x), so that "Hbs can be characterized as the solution space to the system of PDEs: 

62(8) f = e,{d)f = e,{d)f = 0. (14) 

As an ]R[9]-module, I-Lbz generated by the fundamental alternating polynomial 

A-Q^{x) = xiX2X3{xl — x\){x\ — xl){x'\ — x\). (15) 

Since P is at least Vr(A3)-symmetric, Tm\x) must be Vr( A3) -invariant and hence can be 
written in unique ways as polynomials of 62(0;), 63(0;), 64(0;). For m = 1, . . . , 6, one can write 
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Lemma 3.1 For A; = 0, 1, 2, 3, the following hold: 

(1) If none of 0,^2^ ^^i^ "^^ zero, then the infinite system ^ is equivalent to the finite 
system f[T^ so that one has 'Hp(k) = "^Aa- 

(2) If a^^^ = but none of 0^2^ \ ^^e^'' zero, then the infinite system (P) is equivalent to 
the finite system (|T^ so that one has 'Hp(k) = 'Hb^- 
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Proof. First, suppose that none of 0^2^ (^t'^ is zero. A part of equations f lTO]) can then be 
inverted to express 62(0;), 63 (x), 64 (x) as polynomials of t^\x), tI^\x), t'^\x), so that system 
f lT3|) is equivalent to T!^\d)f = T!^^''\d)f = T'^\d)f = 0. For any m > 5 equation Tm\d)f = 
is redundant because Tm\x) is a polynomial of 62 (x), 63 (x), 64 (x). This implies that system 
is equivalent to system (fT5]l . leading to the conclusion of assertion (1). Next, suppose that 
Og^^ = but none of a'^\ a!i \ a^^ is zero. Another part of equations flTBl) can then be inverted 
to express 62 (x), 64 (x), eQ{x) as polynomials of T2''\x), t^''\x), rg'^''(x), so that system (HM is 
equivalent to T!^\d)f = T^\d)f = Tlf^\d)f = 0. Under system ([H]), equation Tm\d)f = is 
redundant for m = 5 or m > 7. Indeed the Vr(A3)-invariance of Tm\x) allows us to write 

2a+3b+4c=m 

with suitable constants a^^),. Equations (HM imply that the summand with index (a, 6, c) 
vanishes if either a > 1, or 6 > 2, or c > 1. Thus the index of a nonzero summand, if any, must 
satisfy a = c = and b < 1 and som = 2a + 36 + 4c<2-0 + 3- l + 4- = 3. Therefore system 
([6]) is equivalent to system f|T^ . This proves the second assertion of the lemma. □ 



4 Skeletons of Triakis Tetrahedra 



Lemma 13.11 is applied to each skeleton of the triakis tetrahedron to establish Theorem 12.11 
4.1 Vertex Problem 

The polynomial Tm\x) in formula ([7j) is divided into two components: 

r^H^) = ri:%{x)+4%{x), (17) 
according to the two types of vertices, where 

Ti%{x) = (A, xr + (B, xr + (c, xr + (d, xr, 

= (a, x)'" + (b, x)"^ + (c, xr + (c, x^. 
Using the coordinates fll2p and the relation d = r(A + B + C)/3 etc., one finds in formulas fll6p . 

4°) = V + 4 > 0, af = ^(3 - r){r^ + 3r + 9), 

af = ^r' + 16 > 0, = ^{r' + 9)(r^ - 9r' + 81) > 0. 

81 24o 

Note that af^ = if and only if r = 3. Thus if r = 3 then assertion (2) of Lemma [3.11 leads to 
the upper case of formula ([5]) with k = 0, while if r 7^ 3 then assertion (1) of that lemma leads 
to the lower case of it. In either case there is no gap between 'Hp(o) and 'Hg(o)- 
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^S(^) = hmiA, AB) + h^{A, AC) + /i™(A, AD) + /i^(B, AB) + h^B, BC) + /i™(B, BD) 
+ hmiC, AC) + /i™(C, BC) + kmiC CD) + /i„(D, AD) + /i^(D, BD) + /i„(D, CD), 

^i,m(a;) = /im(A, Ab) + hm{A, Ac) + hm{A, Ad) + /i^(B, Ba) + /i„(B, Be) + /i„(B, Bd) 
+ hmiC, Ca) + hmiC, Cb) + /i„(C, Cd) + hmiB, Da) + /i„(D, Db) + /i„(D, Dc), 

''"3,m(3^) = ^m(a, Ba) + /im(a, Ca) + hm{a., Da) + /im(b, Ab) + /im(b, Cb) + hm{h, Db) 
+ /i„(c, Ac) + h^ic Be) + hmic, Dc) + /;,™(d. Ad) + /i„(d, Bd) + /i^(d, Cd). 

Table 1: The polynomials ri,m(x), u = 1,2,3, in formula (fT9|) . 
4.2 Edge Problem 

It is obvious that for an index AB G /i, the foot on the line £ab is AB = (A + B)/2, with this 
rule applying to every index of the same type. For indices of the other type in Ji, one finds 

Ah = {(3,-/3, -a), Ac = {/3,-a,-(3), Ad = (a, -(3, -(3), 

Ba= /3,-a), Be = (-a, /3, -/3), Bd = a, -/3), 

Ca = /3), Cb = (-a, /3), Cd = a). 

Da = (a, (3, (3), Dh = {/3, a, /3), Dc = {(3, (3, a). 



where 



4r(r-3) ^ 2r(r + 3) 



3(r2-2r + 9)' 3(r2-2r + 9)' 

For example, formula Ad = {a, —(3, — /3) follows from the condition that the point Ad should 
lie on the line ^Ad, while the vectors Ad and d — A should be orthogonal. 

The polynomial "* {x) in formula (|8]) can be divided into three components: 

i.^) = vei ■ ri^lix) + ve2 ■ r^'^lix) + veg ■ t.^I{x), (19) 

according to the three types (ITU]) of vertex-to-edge incidences, where ri^mix) is given as in Table 
[Hand the abbreviation hm{P,Q) '■= hm{(P,x), (Q,x)) is used for two vectors P, Q G M^. 
The three types of vertex-to-edge incidence numbers are evaluated as 

/ — 9 — T t(t — 1) 

vei = V2, vea = veg = , ^ (20) 

V3(r2-2r + 9) ^3(^2 - 2r + 9) 

Indeed the formula for vei is easy to see. To derive those for ve2 and vea, take a look at the 
edge Ad in Figure O Observing that the unit normal vectors nA,Ad and nd,Ad are given by 
nA,Ad = — 'Md,Ad = (A — d)/|A — d|, where | ■ | denotes the norm of a vector, one can calculate 
ve2 = [A : Ad] = ( A — Ad) / nA,Ad and ves = [d : Ad] = (d — Ad) / nd,Ad as indicated above. 
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Putting all these informations together into formulas f lT6|) . one finds 



a^^^ = 20^2 + ^(r2 + r + 9) v/3(^^^^^27T9)> 0, 
9 



= 96V2 + -(3-r)(r2 + 4r + 9)v/3(7^^27T9), 
9 

= 48^2 + — (r^ + - Sr^ + 9r + 81) v'3(r2 - 2r + 9) > 0, 
81 

64 

= 768^2 + (r^ + - 3r^ - 18r^ - Tlr'^ + 81r + 729) v/3(r2 - 2r + 9). 

243 

Observe that 02"*^^ and positive for every r > 0. Indeed the former is obvious and the 

latter follows from the fact that 'i\}(r) := + — 'ir^ + 9r + 81 has a positive value '?/'(0) = 81 
at r = as well as a positive derivative i\}'{t^ = 4r^ + 3(r — 1)^ + 6 for every r > 0. On the 
other hand, a^r^^ is positive in < r < 3, strictly decreasing in r > 3 and tending to —00 as 
r — i- +00. Thus there exists a unique positive number r = ri at which a^^'^ = 0. Observe that 

,(1) J ^('Q ^^/'^2 I I n^. /q7Z2 I nfi,/o I — ^ (^^2 _ _|_ 3^ Xi(^) 



1^(3 - r){r^ + 4r + 9) v/3(r2 - 2r + 9) _ qq^^ 



so ri must be a positive root of Xil^)) where Xi(^) is the sextic polynomial in Conversely 
one can show that Xil'") certainly has a unique positive root that yields ri = 3.62398 ■ ■ ■ . Note 
that Og^'* = 1661.36 ■ ■ ■ is nonzero at r = ri. Thus if r = ri then assertion (2) of Lemma [3.11 
leads to the upper case of formula with k = 1, while if r 7^ ri then assertion (1) of that 
lemma leads to the lower case of it. There is a gap between 'Hp{i) and T^gci) only when r = ri. 

4.3 Face Problem 

For each index of I2, the foot on the corresponding affine plane is given by 

ABc= (-5,-5,-7), ABd = (5, 5,-7), ACb = (-5,-7,-5), ACd = (5, -7, 5), 
ADb = (7, 5,-5), ADc = (7,-5, 5), BCa = (-7, -5, -5), BCd = (-7, 5, 5), 
BDa= (5, 7,-5), BDc = (-5, 7, 5), CDa = (5, -5, 7), CDb = (-5, 5, 7), 

where 

2r^ ^ r(r - 3) 

7 := 777-^5 ^ — T^' ^ 



3(r2-2r + 3)' ' 3(r2-2r + 3)' 

For example, formula ABd = (5, 5, —7) follows from the condition that the point ABd should 
lie on the plane i^ABd, while the vector ABd should be orthogonal to both d — A and d — B. 
Observe that up to symmetries there are three types of vertex-edge-face flags: 

(1) A ^ AB ^ ABd (VEl) & (EFl), (2) A ^ Ad ^ ABd (VE2) & (EF2), 

(3) d ^ Ad ^ ABd (VE3) & (EF2), 

according to which the polynomial (x) in formula can be divided into three components: 

'^mi^) = vei ■ efi ■ 4!i(x) + ve2 ■ efa ■ 4%{x) + veg ■ efz ■ 4%{x), (21) 
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Ti%{x) = h„,{A, AB, ABc) + /i„(A, AB, ABd) + /i^(A, AC, ACb) + /i„(A, AC, ACd) 
+ /i„(A, AD, ADb) + hm{A, AD, ADc) + /i„(B, AB, ABc) + h^(B, AB, ABd) 
+ h^{B, BC, BCa) + h^{B, BC, BCd) + /i„(B, BD, BDa) + h^{B, BD, BDc) 
+ hm{C, AC, ACb) + h„,{C, AC, ACd) + hm{C, BC, BCa) + /i„(C, BC, BCd) 
+ h^{C, CD, CDa) + /i„(C, CD, CDb) + /i„(D, AD, ADb) + h^{D, AD, ADc) 
+ h^{D, BD, BDa) + /i^(D, BD, BDc) + h^{D, CD, CDa) + /i„(D, CD, CDb), 

'^2%i^) = ^m(A, Ab, ACb) + /im(A, Ab, ADb) + /im(A, Ac, ABc) + /im(A, Ac, ADc) 
+ h^{A, Ad, ABd) + /i„(A, Ad, ACd) + /i™(B, Ba, BCa) + hmiB, Ba, BDa) 
+ hmiB, Be, ABc) + hmiB, Be, BDc) + /i„(B, Bd, ABd) + /i„(B, Bd, BCd) 
+ hmiC, Ca, BCa) + h^iC, Ca, CDa) + /i„(C, Cb, ACb) + /i„(C, Cb, CDb) 
+ hmiC, Cd, ACd) + hmiC, Cd, BCd) + /i^(D, Da, BDa) + /i„(D, Da, CDa) 
+ /i^(D, Db, ADb) + hmiB, Db, CDb) + /i„(D, Dc, ADc) + /i^(D, Dc, BDc), 

r3^2i(^) = hmia., Ba, BCa) + hmia., Ba, BDa) + /imla, Ca, BCa) + Ca, CDa) 
+ hmia, Da, BDa) + hmia. Da, CDa) + /^^(b, Ab, ACb) + /^^(b, Ab, ADb) 
+ hmih, Cb, ACb) + /i„(b, Cb, CDb) + /i„(b, Db, ADb) + /i^(b, Db, CDb) 
+ /im(c, Ac, ABc) + hmic, Ac, ADc) + hmic. Be, ABc) + hmic. Be, BDc) 
+ hmic, Dc, ADc) + /i„(c, Dc, BDc) + /i^(d. Ad, ABd) + /i„(d. Ad, ACd) 
+ hmid, Bd, ABd) + hmid, Bd, BCd) + /i^(d, Cd, ACd) + /i„(d, Cd, BCd). 

Table 2: The polynomials Ti^^mix), v = 1,2,3, in formula (12T|) . 

where the polynomials rifm(x), = 1,2,3, are given as in Table [2] and the abbreviat'^^ 
hmiP,Q,,^) '■= hmi{P,x), {Q,,x), (R, a;)) is used for three vectors P, Q, R e M^. 

While the vertex-to-edge incidence numbers are given in fl2Ul) . the edge-to-face ones are 

3 - r . v^r(r - 1) 

efi = , , „ . , ef2 - 



v/3(r2-2r + 3)' v^(r2 - 2r + 3)(r2 - 2r + 9) ' 

So upon multiplying by a nonzero constant simultaneously, one may put 

„ 3 — r (9 — r)r(r — 1) 

vei ■ efi = — — - — — , vea ■ efa - 



r2-2r + 3' (r2 - 2r + 3)(r2 - 2r + 9) ' 

r2fr — 1)2 



vea ■ ef 2 



(r2-2r + 3)(r2-2r + 9)' 



Notice that what is important in expression ( 12T]) is only the ratio vei ■ efi : ve2 ■ ef2 : ves ■ ef 
Putting all these informations together into formulas ( IT6|) . one finds 

af) = ^(r^ + 2r + 15) > 0, af^ = y (3 - r)ir^ + 5r + 12), 

32 128 

af ) = (r^ + 2r3 - 3r^ + 81) > 0, = (r*^ + 2r^ - 3r^ - 27r^ - 5Ar^ + 81r + 972) 

27 81 
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c+ 




Figure 4: Adjoining a pyramid to a face of a regular octahedron; OD'*'"^ : Od"'"^ = r : 1. 

(2) (2) (2) 

It is easy to see that 02 and 04 are positive for every r G M and that 03 = if and only 
if r = 3, in which case ag^"* = 1536 is nonzero. Thus if r = 3 then assertion (2) of Lemma 13.11 
yields the upper case of formula ([5]) with k = 2, while if r 7^ 3 then assertion (1) of that lemma 
yields the lower case of it. In either case there is no gap between 'Hp(2) and 'Hg(2)- 

By O Theorem 2.2] the volume problem {k = 3) has the same solution as the face problem 
(A; = 2), since P is isohedral. The proof of Theorem 12. II is now complete. 

5 A Family of Isohedral Triakis Octahedra 

Let f2 be a regular octahedron with center at the origin O. The six vertices of Q can be written 
A^, B"^, C^, where A"*" and A~ are antipodal to each other and so on. The eight faces are 
then given by A'^B^C^ with a,b,c = ±. An isohedral triakis octahedron P is obtained from Q 
by adjoining to each face of it a pyramid of appropriate height, or excavating such a pyramid, 
just as in the construction of an isohedral triakis tetrahedron. Let D"'"^ denote the top vertex 
of the pyramid based on the face A"B^C^. The polyhedron P depends on a parameter r > 
in such a manner that the distance ratio OD'^''^ : Od"^'^ is r : 1, where d"'"^ is the center of the 
face A'^B^'C^ (see Figure S]). 

Let us look more closely at the polyhedron P for various values of r. The values r = 1 
and r = 3/2 are special in that as a point set, P degenerates to the original octahedron Q at 
r = 1 and to a rhombic dodecahedron at r = 3/2 where two neighboring faces, say, A+B+D+++ 
and A+B+D^"*"" are coplanar. Note that P is convex if and only if 1 < r < 3/2, in which 
interval the value r = 3(-\/2 — 1) = 1.24264- ■ ■ is distinguished in that P becomes a Catalan 
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Figure 5: The isohedral triakis octahedron with r = 1.82977 ■ • • . 

sohd [HIZ], or an Archimedean dual sohd, called the small triakis octahedron, whose dual is the 
truncated cube [3, 8, 8]. As in the case of triakis tetrahedron we employ the convention that as 
a combinatorial polyhedron, P has the constant skeletal structure for all r > 0, even at r = 1 
and r = 3/2. Figure O exhibits the shape of P when r = 1.82977 ■ ■ • ; see the k = 2,3 case of 
Theorem 15. II for the origin of this particular value. 
If the vertices of our octahedron Q are taken as 

A± = (±1,0,0), B± = (0,±1,0), C± = (0,0,±1), (22) 

then the symmetry group of Q is the same as that of the cube C in ^ namely, the Weyl group 
1^(63). It is just the symmetry group G = G{k) of P{k) for every A; = 0, 1, 2, 3 and r > 0. 

Theorem 5.1 For the family of isohedral triakis octahedra, 'HG{k) £ 'Hp(fc) if o,nd only if 

• k = andr = ro:=3- 2~^l^ = 1.78381 ■ ■ ■ , 

• k = 1 and r = ri := 2.24580 ■ ■ ■ is the unique positive root of an octic equation 

Xi(r) := 16r^ + 32r^ + 40r^ - 48r^ - 396r^ - 432r^ - 810r^ - 972r - 729 = 0, (23) 

• k = 2,3 and r = r2 = := 1.82977 ■ ■ ■ is the unique positive root of a quartic equation 

X2{r) := 4r^ + 8r^ + 6r^ - 18r - 81 = 0. (24) 

// any of these is the case, then dim'HG{k) = 48 < dim'Hp(^k) = 96 and as W[d]-modules, 
'HG{k) = "Hbs "is generated by the polynomial Z\B3(a;) in formula fll5p while 'Hp(k) is generated 
by 

F{x) := Z\B3(a;) ■ {5(x^ + 4 + x^) - 13(a;^a;^ + xjxl + xlxl)}. (25) 

As in the case of Theorem 12 . 1 1 this theorem is also established through the analysis of PDEs 
(E]). To adjust the notation to the current situation we represent the index sets Jq, /i, I2 in §3] 
by letting the vertices, edges and faces of P to speak of themselves, that is, by setting 

/o = {A", B^ C"} U {D"^^}, 

Ji = {A"B^ A"C^ B''C^} U {A"D"''^ B^D"^^ C^D"^"}, 
I2 = {A'^B^D'^^^ A^C'D'^^^ B^'C'D"''^}, 
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where a, b, c take ± signs freely and {■ ■ ■ } stands for an orbit under symmetries. For an index 
A"B* G /i the same symbol A"B^ denotes the foot of orthogonal projection from the origin O 
to the affine line passing through A" and B''; this rule also applies to another index A^D'^''^ G /i 
as well as to all the other indices of Ji. Similarly, for an index A^B^D'*'"^ G I2 the same symbol 
y^agfe-Qafec (jgnotes the foot of orthogonal projection from O to the affine plane passing through 
A", B^, D"'"^, with this rule applying to all the other indices of I2', see Figure IH 
Up to symmetries there are three types of vertex-to-edge incidence relations: 

(VEl) A+ ^ A+B+, (VE2) A+ ^ A+D+++, (VE3) D+++ ^ A+D+++. 

On the other hand, up to symmetries there are two types of edge-to-face incidence relations: 

(EFl) A+B+ -< A+B+D+++, (EF2) A+D+++ ^ A+B+D+++. 

The incidence number of an incidence relation depends only on its type. Let ve^, denote the 
incidence number of type (VEz/), u = 1,2,3, and efi, denote that of type (EFu), u = 1,2, 
respectively. Notice that the presentation around here is quite similar to the presentation 
around formulas (ITU]) and (ITT]) . Indeed the treatment of triakis octahedra will be very parallel 
to that of triakis tetrahedra, so that the main focus in what follows will be on how one can 
modify the arguments in the tetrahedral case to the octahedral one. 

We consider how Lemma 13.11 can be modified. Since P is iy(B3)-symmetric, must 
be iy(B3)-invariant and hence can be written in unique ways as polynomials of 62(0;), e4(x), 
eQ{x). Note that Tm\x) is identically zero whenever m is odd. For m = 2, 4, 6, 8, one can write 
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Lemma 5.2 For = 0, 1, 2, 3, the following hold: 

(1) If none of 0^2^ \ (^6^^ ^-^ zero, then the infinite system ([6]) is equivalent to the finite 
system ( 1T4|) so that one has 'Hp(k) = "Hbs- 

(2) If a^^^ = but none of a'^\ a^Q \ a^p is zero, then system ([6]) is equivalent to 

e2{d)f = e,{d)f = el{d)f = Q. (27) 
If S denotes the solution space to system ( 127|) . then there exists an exact sequence 

, s ^ nB, — > 0. (28) 

As an M.[d]-module, S is generated by the polynomial F{x) defined in formula ( 125|) . 
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Proof. The proof of assertion (1) of this lemma is the same as that of assertion (2) of Lemma 
13.11 Under the assumption of assertion (2) of this lemma, the first, third and fourth equations 
of (l26l) imply that system (127|) leads to T!^\d)f = TQ''\d)f = Tg''\d)f = and conversely the 
latter leads back to the former. Under system fl27|) . equation Tm\d)f = is redundant for 

(k) 

m = 4 or for any even m > 10. Indeed, for m = 4 it is trivial from a\ = and the second 
equation of fl26l) . For m > 10 the 14^(B3)-invariance of Tm\x) allows us to write 

r^\d)f= J2 ^a^c-e2id)-el{d)-el{d)f, 

2a+4b+6c=m 

with suitable constants a^^).. Equations (1271) imply that the summand with index (a, b, c) 
vanishes if either a > 1, or 6 > 2, or c > 1. Thus the index of a nonzero summand, if any, must 
satisfy a = c = and b < 1 and som = 2a + 46 + 6c<2-0 + 4- l + 4- = 4. Therefore 
system is equivalent to system (1271) . Next we verify the exact sequence (1251) . The inclusion 
"Hbs C S and the well-definedness of 64(9) : S — )• "Hbs are obvious since "Hbs is the solution 
space to system ([Hj) while S is to system (127|) . For the same reason sequence (l28l) is exact at 
the middle term S. A direct check shows that polynomial F{x) in formula (125|) satisfies 

e2{d)F = ee{d)F = 0, e4(9)F = -15120Z^B3- 

Hence F E S and 64(9) sends F to Ab^ up to a nonzero constant multiple. Thus 64(8) : S — )■ 
"Hbs is surjective, since it is an R[(9]-homomorphism and Ab^ generates "Hbs as an M[9]-module. 
Finally we show that F generates S as an M[(9]-module. For any / G 5, consider 64(9)/ G "Hbs- 
There is a polynomial such that 64(8)/ = —15120 v9(9)Z\b3- Then g '■= f — f{d)F G 5 
satisfies e4{d)g = 64(8)/ - ^{d)e4{d)F = 64(8)/ + 15120 (^(9)Z\b3 = 0. Exact sequence 
tells us that g G T-Lb^ and so there is a polynomial ip{x) such that g = — 15120 '0(5)Z\b3 = 
ilj{d)e4{d)F. Now if ri{x) := V5(x) + ilj{x)e4{x) then / = ri{d)F. This proves the last claim. □ 



6 Skeletons of Triakis Octahedra 

Lemma 15.21 is applied to each skeleton of the triakis octahedron to establish Theorem 15.1 



6.1 Vertex Problem 

Formula (1T7|) for triakis tetrahedra carries over to the case of triakis octahedra, but this time 
formula (ITS]) should be replaced by 



5^(A^ xr + + ^r, A%{x) = 5^ {d^"^, xy 

a b c {a.,b,c) 



where the sums are taken over all a,b,c = ±. Using the coordinates ( 12 2 p and the relation 
jyabc ^ ^(^^a + + C'=)/3, one finds in formulas 



(0) 8 2 , o ^ n (0) 32 / 4 3 

On ' = -r + 2 > 0, a\' = — r 

^ 9 ' ^ 81 V 23 

af = ^{Ar^ + 9)(16r^ - 2>Qr^ + 81) > 0, = + 4 > 0. 

243 d5o1 
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Note that a^^^ = if and only if r = tq := 3 ■ 2^'^/^ = 1.78381 ■ ■ ■ . Thus if r = tq then assertion 
(2) of Lemma 15.21 yields 'Hp(o) = S, while if r 7^ tq then assertion (1) of Lemma 15.21 yields 
'Hp(o) = "Hbs- One has "Hgco) ^ 'Hp(o) only when r = tq. This proves Theorem 15.11 for k = 0. 

6.2 Edge Problem 

It is obvious that for an index A°B^ G Ji, the corresponding foot is A^B^ = (A'' + B*)/2 with 
this rule applying to every index of the same type. For indices of the other type in Ji, one finds 

j^aj^abc ^ (^^^ ^^^^ ^bj^abc ^ ^^^^ ^ ^cj^abc ^ ^^^^ ^ 

with a, 6, c = ±, where 

2r^ r(3 - r) 

u := ——z -, V ■— 



3(r2-2r + 3)' ' 3(r2-2r + 3)' 

Formula f lT^ remains true for triakis octahedra if formulas in Table [T] are replaced by 

+ J2 ^-(B', B'C^) + J2 ^-(C^ A°C^) + 5^ hUC', B'C'), 
r^li^) = Yl ^-(A''' + E ^™(B^ B^'D'^^^) + J2 hmiC', CT"''^), 

The three types of vertex-to-edge incidence numbers are evaluated as 

1 3 — r r(r — 1) 

vei = ve2 = , , vea = , (29) 

V2 ^3(r2_2r + 3) y/3{r^ - 2r + 3) 

Putting all these informations together into formulas ( 12^]) . one finds 



a 



(1) 
2 
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V2 + -{r^ + r + 3) V3(r2 _ 2r + 3) > 0, 



9 



16 

= -12^2 + — (2r^ + 2r^-9r-27)v/3p^^27T3), 
81 

= 12^2 + — (16r^ + 16r^ -18r^ + 81r + 243) v/3p^^27T3)>0, 



a 



(1) 



16 

I2V2 + (8r^ + 8r^ - 36r^ - 108r^ - 162r3 + 729r + 2187) v/3(r2 - 2r + 3). 



6561 



Observe that 03 '^e ^'^^ positive for every r > 0. Indeed the former is obvious and the 

latter follows from the fact that '^(r) := 16r^ + 16r^ — 18r'^ + 81r + 243 has a positive value 
ip{0) = 243 at r = and a positive derivative ip'{r) = 96r^ + 71r^ + {Sr"^ — 9)^ for every r > 0. 
On the other hand, there exists a unique positive number r = ri at which 04^'' = 0, because 



a 



da^P imr^ir"^ - r + 11 



■4 



-4(4 + 3^2) < at r = 0; — ^ = , ^ L > for r > 0, 

dr 27v/3(r2 - 2r + 3) 
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and a^^^ tends to +00 as r — )■ +00. If Xi(^) is the octic polynomial defined in formula (12^ . 



|l2V2 + i^(2r^ + - 9r - 27) v/3(r2 - 2r + 3)| 



2187 



(2r^-4r + 3) 



(1) 

54.1247- ■ • is nonzero at r = ri. Thus Lemma [5.21 leads to Theorem 15. II for k = 1. 



So ri = 2.24580 ■ ■ ■ is the unique positive root of octic equation Xii^") = 0- Note that ag 



6.3 Face Problem 

For each index of I2, the foot on the corresponding affine plane is given by 

j^a^b^abc ^ ^^^^ ^^^^ pa(^c^abc ^ ^^^^ ^ B^C'^D'^''" = {ap, bq, cq) , 

with a, 6, c = ±, where 



r(3 — 2r) r 



2 



P--=^T77r-o T— Q--- 



3(2r2-4r + 3)' ^' 3(2r2-4r + 3)' 
Observe that up to symmetries there are three types of vertex-edge-face flags: 

(1) A+ ^ A+B+ ^ A+B+D+++ (VEl) & (EFl), 

(2) A+ -< A+D+++ ^ A+B+D+++ (VE2) & (EF2), 

(3) D+++ ^ A+D+++ ^ A+B+D+++ (VE3) & (EF2). 
Formula fl2T]) remains true for triakis octahedra if formulas in Table [2] are replaced by 

-(2) _ J2 /^™(A^ A''B^ A^B^'D'^''^) + ■ ■ ■ sum over flags of type (1), 



J2 hra{^-\ A'^D'^^^ A'^B^D"'"^) + • • • sum over flags of type (2), 



^S(^) = ^m(D'''", A'^D"^^ A"B''D"''^) + ■ ■ ■ sum over flags of type (3). 

While the vertex-to-edge incidence numbers are given in fl29|) . the edge-to-face ones are 



3 - 2r . r{r - I] 

eii = , , ^ ef2 



v/6(2r2-4r + 3)' v^(r2 - 2r + 3)(2r2 - 4r + 3) ' 

So upon multiplying by a nonzero constant simultaneously, one may put 

„ 3 — 2r (3 — r)r(r — 1) 

vei ■ ef 1 = — —7, ve2 ■ ef2 - 



2(2r2-4r + 3)' (r2 - 2r + 3)(2r2 - 4r + 3) ' 

r2(r — 1)2 
(r2-2r + 3)(2r2-4r + 3)' 

Putting all these informations together into formulas f l26p . one finds 
= ^(r^ + 2r + 6) > 0, af = ^ X2(r), 



= —{IQr^ + 32r^ + 24r^ - 18r^ - 54r^ + 243) > 0, 
81 

,(2) 
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-(16r^ + 32r^ + 24r6 - 72r^ - 324r^ - 648r3 - 486r2 + 1458r + 6561), 
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where X2(^) is the quartic polynomial defined in (^^. Observe that 05 and are positive 

for every r > 0. Indeed the former is obvious and latter follows from the fact that a^^ {r > 0) 

(2) 

attains its minimum 22.0304 ■ ■ ■ > at r = 0.743471 ■ ■ ■ . Note that 04 = if and only if r 
is the unique positive root r2 = 1.82977- ■■ of equation (124|) . at which = 13.2853 ■• ■ is 
nonzero. Lemma [5.21 thus leads to Theorem 15.11 for k = 2. 

By |S1 Theorem 2.2] the volume problem {k = 3) has the same solution as the face problem 
(A; = 2), since P is isohedral. The proof of Theorem 15.11 is now complete. 

It is an interesting exercise left behind to deal with a similar problem for isohedral triakis 
icosahedra. It would also be interesting if such special solids as discussed in this article appear 
in nature and physical sciences. 
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